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Abstract 

In this Letter we discuss theoreticaUy the instabihties of thermotropic freely 

standing smectic-A films under shear flow We show that, in Couette 

geometry, the centrifugal force pushes the liquid crystal toward the outer 

boundary and induces smectic layer dilation close to the outer boundary. 

Under strong shear, this effect induces a layer buckling instability. The critical 

shear rate is proportional to \/^fd, where d is the thickness of the film. 
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Consider shear flow experiments on freely standing fluid films in a Couette geometry 
(Fig. 1), in which the inner cylinder is rotating with speed Vq, and the outer cylinder is at rest. 
Both cylinders are perpendicular to the plane of the film. When a soap film is used in the 
experiment, the laminar flow between the cylinders is stable up to very large vq, and the film 
gets thicker towards the outer cylinder due to centrifugal force 0]. Brightly colored patterns 
are exhibited because of this thickness variation. However, for smectic-A films with smectic 
layers parallel to the air-fluid interface, which were treated as two-dimensional systems in 
many analyses [§], there exists a critical speed Vc, such that when vq is above Vc, defects are 
generated This clearly indicates that the layer structure becomes important, and the 
system can no longer be treated as two-dimensional. On the other hand, experimental 0-0] 
and theoretical studies on the rheological properties of bulk smectic-A systems show 
that under shear flow, bulk systems are most stable when the layers are parallel to the shear 
plane. Thus the finite thickness and/or the free surfaces have to be essential to the observed 
defect generation in the experiments on freely standing smectic-A films. 

In this Letter we study the rheology of freely standing thermotropic smectic-A films 
under shear flow. In a Couette geometry the centrifugal force pushes the liquid crystal 
away from the inner wall, which induces a layer dilation close to the outer wall. When 
this dilation exceeds a certain limit, the smectic layers are no longer stable against small 
undulational perturbations, hence defects are generated. This instability is similar to the 
classical Helfrich-Surault type instabilities but the film has free surfaces and is far from 
equilibrium. For typical materials and experimental setups the strain induced by centrifugal 
force provides the critical velocity, Vc, comparable to the experimentally observed result 
Thus we conclude that this scenario is a good candidate for the mechanism which is 
responsible for the defect generation in Ref. [Q]. Our analysis indicates that this instability 
occurs in freely standing smectic-A films but not in bulk smectic-A, since it is a result of 
the interplay between crossover and boundary effects (2D to 3D and the existence of free 
surfaces), and the external driving flow. 
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Strain induced buckling instability for freely standing smectic-A films. — We begin with 
the classical strain induced buckling instability for smectic-A in a new geometry, i.e., freely- 
standing films. Let the smectic layers be parallel to the surfaces of the film, which are 
located ai z = ±d/2, and ignore the boundaries in the x and y directions. The total elastic 
free energy of the film is [ pT| , P^ 



F = a(^J dA+ + j dA^ 

+-/ dh^ j ^^dz(B{E[u]f + K,{dluf) , (1) 



where / dA± = J d'^r±yl + (d±[u]z=±d/2j is the area of the surfaces, u is the smectic layer 
displacement, a is the surface tension, Ki and B are layer undulation and compression 
elastic moduli, and E[u] = d^u — , is the layer dilation. The continuity of the 

normal stress at the surfaces requires ||T3|-[T5| 



± [BE[u]l^^,^, = a{6A^/6u) + P.^, (2) 

where S/5u stands for variational differentiation, are the external normal stresses acting 
on the surface at z = ±d/2. 

Suppose the smectic layers are subject to a uniform but opposite external stress = 
±Pext on the surfaces; for small Pext the layers are uniformly dilated with d^u = Pext/B = e . 
We study the stability of small layer undulations of the form u = uq + ui with uq = ez and 
Ml = UN{qy) sm{qyy) cos{qzz), where, for given g^, > 0, > is chosen such that Eq. (||) is 



satisfied, and the normalization factor N{qy) is chosen such that [N{qy)]^ I-d/2 cos'^{q^ 



d/2 ' 



1 . 

Expanding F to 0{U'^) and using the boundary condition for ui to the corresponding 
order, one finds 

1 r. 
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B{ql - eq') + K^q^ U\ (3) 



where, from Eq. (0), q^ satisfies the following equation. 



Qzd = {Idqyf cot ^ , (4) 



in which Id = J ad/ B is a. characteristic length. As discussed in Ref. |T^, this length char- 



acterizes the crossover from behavior dominated by surface tension to behavior dominated 
by layer elasticity in freely standing smectic-A films. Since we are interested in the onset of 
instability, we only consider the smallest Qz which satisfies Eq. (§) throughout this Letter. 
From Eq.(D it is clear that when the induced strain e is large enough, there is an instability, 
and one would expect growth of the amplitude U. 



From Eq.(^), the critical strain for given Qy is e^Qy) = Qy + ^'^Qy , where X = J Ki/ B is 



a characteristic length for smectic-A ||TT[. Defining dimensionless parameters X = q^d , Q = 
qyld , then 



where g = aj \JKxB is a dimensionless parameter; typically g ~ 0(1). The exact value of 
X cannot be found analytically, but numerically one finds that T{Q,g) has a minimum at 
Q = Q* > 0. From the definition of T{Q,g), Q* depends only on g. For given material 
parameters, one can find the minimum critical strain , e* = ec(Q*), numerically, which 
depends only on g and \/d. Fig. 2 shows Q* and 6e*d/Xg as functions of g. Notice that 
there is a gm such that when g < gm, the first unstable mode has infinite wavelength. 
We can estimate this parameter by expanding the equation for X around Q = 0, and 
then substituting this approximate value of X into Eq. (|^). We find X^/Q"^ + Q^/g^ = 
2 + (1/^2 _ i/3)g2 ^ 0{Q^). For the first unstable mode to occur at Q = 0, the coefficient 
of the term in the small Q expansion of the critical strain has to be positive, yielding 
dm = V^- On the other hand, when 1/g = \JKxB ja — » 0, the surface tension is sufficiently 
large such that the instability essentially becomes the Helfrich-Surault instability originally 



derived for hard boundaries, where the first unstable mode is g* = y vr/dA and the minimum 
critical strain is 2'K\/d. |T0[ 

In the shear flow experiments conducted in Ref. the freely standing smectic-A films 
were made of 4 — cyano — 4' — octylbiphenyl (8CB), which has a ~ 25 dyn/cm |]TB[, B ^ 



5 X 10" dyn/cm |T^, fCi ^ 5 x 10 "dyn [T^, i.e., g ^ 5. From Fig. 2 one finds that (in 



equilibrium), under a uniform layer dilation, the minimum critical strain e* ^ Xg / d = a/Bd, 
and the first unstable mode has wavelength on the order of la- 
Shear induced instability for freely- standing smectic-A films in the Couette geometry. — 
The experimental setup is shown in Fig. 1. We choose cylindrical polar coordinates to 
describe the system in which the z axis is the axis of the cylinders. The radii of the inner 
and outer cylinders are i?/, and i?o, respectively. Let the film be suspended horizontally 
on the rip plane between the inner and outer cylinders with thickness d (from z = —d/2 
to z = d/2). |]l9l In principle d depends on position in the r0 plane. The smectic layers 
are parallel to the r0 plane away from the meniscus. The films are usually thicker near the 



meniscus due to wetting pO|], and the liquid crystal near the meniscus acts like reservoir 



of material in shear flow experiments. Furthermore, the interactions between the liquid 



crystal and the walls also cause edge dislocations ||2T|. Since our focus is the effect of shear 
flow, to avoid additional complications, we assume both walls to be neutral to the liquid 
crystals, and neglect the effect of gravity and the surrounding air so that in equilibrium the 
smectic film is perfectly fiat with a uniform thickness d = d^. As a result, in our model there 
is no material reservoir near the meniscus, and the total volume of the film under shear flow 
is the same as in the film at rest, a simplification which does not occur in real experiments. 



The equations of motion for this system in the isothermal, incompressible limit are [|T2 



du ^ , , 

— + V ■ Vm = f2 + Cp/i , (6) 
p + V ■ Vv^ = -Vp + /iz + r/VV , (7) 

where v also satisfies the incompressibility condition V ■ v = 0. Here v is the velocity of 
the hquid crystal, Cp is the permeation constant, p is the pressure, and h = d ■ {6F/6du). In 
general the viscosity rj should be a second rank tensor due to the anisotropy of the smectic- 



A phase [^,0. However, to simplify the algebra we choose ?7 to be a scalar. Since the 
imposed flow direction is parallel to the smectic layers, and the system is isotropic in the 
plane defined by the layers, this simplification should not affect the qualitative result of our 
calculation. 



The stress should be continuous across the free surfaces, which leads to the following 



conditions, 



l-p + BE[u]Tadlu + 2r]^\ =0. (10) 

V J z=±d/2 

We take no-slip boundary conditions for the velocity on the inner and outer walls, i.e., 
V = t>o at r = i?/, and v = at r = Rq. We first consider the steady state. Subsequently 
we will study the linear stability of the steady state. We assume rotational invariance around 
the z axis. As a result it is natural to set Vr = everywhere. Then incompressibility and 
no-slip boundary conditions lead to = 0. The momentum equation, Eq.(|^, now becomes 

(U) 

r par 

r dr dr dz^ 
0--| + ft. (13) 

The solution of Eq.(|l2|) which satisfies the boundary conditions is f ^ = ar + ^, where 

With Vz = and v^f, given above, the steady state layer displacement has to satisfy (from 
Eqs.dl) and (|13D) h = Bd^u — Kid\u = 0, and ff = 0, i.e., the pressure depends only on r. 
Now one can solve Eq.(pJ]) by direct integration, leading to 

p{r) -p{Ri) 

^,(^(.= -fl?).2.Mn^4(i,-i)) . (14) 

where p{Ri) will be determined later, by conservation of total volume. Since p{r) — p{Ri) 
is positive definite for non- vanishing f^, the pressure near the outer wall is larger than the 
pressure near the inner wall. The boundary condition for the normal stress, the condition 
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h = 0, and Eq.([l4D together with the conservation of total volume will completely determine 
the steady state configuration. 

To understand the effect of this centrifugal force induced pressure, let us first consider 
the long wavelength behavior of the layer displacement to linear order in u. To this order 
the boundary condition for the normal stress is B [dzu]z=±d/2 ~ l^i.'^]z=±d/2 + [p\z=±d/2 ■ 
Simple dimensional analysis (replacing dz by 1/d, d± by l/{Ro — Ri), and typically B/d^ 
a/{Ro — RiY EH) shows that, for nonvanishing p, the surface tension term is negligible 
compared to the layer compression for typical experiments, and the layer dilation on the 
surfaces is related to the pressure by 

B [dzu]^^^^/^ = [p],=±rf/2 • (15) 

The pressure "pushes" the smectic layers on the free surfaces, but with a very large char- 
acteristic in-plane length scale {Rq — Ri)- Since the walls are neutral to the liquid crystal, 
and h = Q everywhere, the layer displacement is approximately u = z = ez . This 
layer displacement causes a change of the smectic film thickness from do to (io(l +p{r)/B). 
Since the total volume of the liquid crystal in the steady state is the same as that of the 
initial state; we determine the pressure at r = Rj from the condition jj^^ dr{27rrdQ) = 

One finds the layers are compressed for small r but are dilated 
for large r due to the effect of centrifugal force and conservation of total volume. 

Consider a small patch located near the outer cylinder having size I ^ Rq — Ri- The 
stability to a layer undulation with wave vector in the radial direction can be analyzed as 
follows. Choosing x = —(f>, y = f, when r is close to Rq with Rq — r ^ /, we can work in 
Cartesian coordinates and ignore the effect of the walls. The small layer undulation can be 
expressed as UN{qy) sm{qyy) cos{qzz). Since the shear flow is in the x direction, the small 
layer undulation with in-plane wavevector in the y is not convected by the flow (i.e., in 
Eq.(^, V ■ Vm vanishes), and it evolves towards h = 'V ■ 5F/5'Wu = 0. Thus the equilibrium 
analysis can be used to show that for large enough vq, the induced layer dilation close to 
the outer cylinder can cause a layer undulation instability. The fact that the flrst unstable 
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JI^° dr hnrd 



Pjr) 
B 



+ 1 



mode for 8CB has a wavelength on the order of Id justifies our use of Rq — Ri ^ I . 

We are now in a position to estimate Vc for the experimental geometry used Ref []I|, i.e., 
Rj = 1.2 cm, Ro = 1.5 cm. Choosing r ~ Rq, and estimating e* ~ a/Bd to estimate 
the layer dilation close to the outer cylinder and the critical strain, one finds, for a typical 
material, p ^ 1 g/cm^, a ^ 25 erg/cm^, f ^ ~ 1500 cm^/sec^. For d = 3.5 x 10~^cm, one 
finds fc ~ 2.1 X lO^cm/sec, and the critical shear rate is 7c = Vc/ {Ro — Ri) ~ 7 x lO'^sec"^. 
The experimental critical shear rate is about 10^ sec~^. Considering the simplifications 
made in this model, the agreement indeed indicates that the effect of layer dilation induced 
by centrifugal force can cause the instability observed in the experiments. Notice that our 
model shows the initial instability should occur close to the outer boundary of the Couette 
cell, and this is consistent with the experimental observation. 



Concluding Remarks. — We have shown that freely standing smetic-A films are in prin- 



ciple unstable against strong shear flow. The characteristic length Id introduced in Ref. |T4 
plays an important role. We considered specifically the Couette geometry and included the 
effects of centrifugal force. We showed that centrifugal force is capable of inducing a layer 
dilation close to the outer cylinder. When the shear rate is large enough, defects are gener- 
ated due to a layer buckling instability similar to Helfrich-Surault type instability |TD|, and 



the calculated critical shear rate is on the same order as the experimental measurements. 

However, notice that the linear stability analysis in this Letter is appropriate only for 
layer undulations with in-plane wavelength small compared to Rq — Ri- As we point out, 
when g = a/y/KiB < the first unstable mode against a uniform layer dilation has very 
long wavelength. Hence, our analysis does not apply to certain materials with weak surface 
tension. For the material (8CB) used in Ref. flj, (7 ~ 5, and our analysis is appropriate. 

Our calculations have assumed perfectly aligned smectic layers and ignored any interac- 
tions between the air or solid boundaries with the liquid crystal. Edge dislocations and the 
material reservoir commonly observed near the meniscus I^T], which may not be negligible 
in the experiment, are also ignored in our approach. Due to these simplifications, our work 
should only be compared semi- quantitatively with the experiments. We also note that, un- 

8 



der strong shear flow , defects are generated in freely standing lyotropic smectic-A films 
in Couette geometry. Since the density of the solvent is an extra hydrodynamic variable in 
lyotropic systems, p^|25| whether the mechanism for defect generation is the same as that 
for the thermotropic systems remains to be answered. However, our study has provided 
the basic ingredients of the experimental instability, i.e., the finite film thickness, surface 
tension, and the geometry of the experimental setup. It also provides the physical picture 
of the instability, i.e., the differences between a smectic-A film and a soap film in Couette 
geometry, and the differences between a bulk system and a freely standing film. Hence we 
believe that we have achieved our goal, which was to elucidate the basic mechanism behind 
the interesting experiments reported in Ref. We also add to the growing catalog of 
fascinating properties of smectic-A films far from equilibrium. 

We thank D. Dash and Professor X-1. Wu for very helpful discussions. Professor T. Ohta 
was very helpful in the early stages of this work, and we thank him. D.J. and H.-Y. C. are 
grateful for the support of the NSF under Constract No. DMR9217935. 
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FIGURE CAPTIONS 



Figure 1. Schematic of the experimental setup in Couette geometry, the z axis points out 
of the paper 

Figure 2. The first unstable mode Q* and the critical strain Qe*d/Xg for freely standing 
smectic-A films under a uniform strain. 
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